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ABSTRACT 


We  present  infrared  (IR)  absorption  spectra  over  the  800-7800  cm'1  region  of  cryogenic 
parahydrogen  (pHi)  solids  doped  with  H2O,  D2O  and  HDO  molecules.  Analysis  of  the 
rovibrational  spectra  of  the  isolated  H2O,  D2O  and  HDO  monomers  reveals  their  existence  as 
very  slightly  hindered  rotors,  typically  showing  only  2  to  5  %  reductions  in  rotational  constants 
relative  to  the  gas  phase.  The  nuclear  spin  conversion  (NSC)  of  metastable  J  =  1  ortho-H20 
(0H2O)  and  para^O  (PD2O)  molecules  follow  first  order  kinetics,  with  single  exponential 
decay  lifetimes  at  T  =  2.4  K  of  1900±100  s,  and  860±50  s,  respectively.  We  report  without 
discussion  some  absorptions  of  water  clusters  produced  during  sample  annealing.  We  report  and 
assign  a  number  of  absorptions  to  oH2-water  pairs  or  “complexes.”  The  main  features  of  the 
0H2-H2O  and  0H2-D2O  spectra  are  explained  qualitatively  by  assuming  a  semi-rigid  C2V  structure 
with  the  0H2  acting  as  a  proton  donor  to  the  O  atom.  Surprisingly,  NSC  of  oH2-water  complexes 
proceeds  at  very  nearly  the  same  rate  as  for  the  corresponding  water  monomer.  We  report 
unassigned  spectra  of  larger  (oH2)n-water  clusters,  and  the  even  more  surprising  observation  of 
the  prolonged  survival  of  oH20  and  pD20  molecules  clustered  with  several  oH2  molecules.  We 
report  and  assign  a  number  of  water  dopant-induced  IR  absorption  features  of  the  pH2  host,  along 
with  cooperative  water-pH2  transitions  in  which  the  vibrational  excitation  of  the  pH2  solid  is 
accompanied  by  a  pure  rotational  transition  of  the  water  dopant. 

Keywords: 

matrix  isolation  spectroscopy;  solid  parahydrogen;  nuclear  spin  conversion;  dopant-induced  IR 
activity;  cooperative  absorption 
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I.  INTRODUCTION 

The  Winnewissers’  solid  hydrogen  spectroscopic  opus  (1-6)  has  yielded  the  discovery  of 
several  new  rovibrational  transitions,  including  the  novel  class  of  condensed  phase  “triple 
transitions,”  in  which  three  hydrogen  molecules  cooperate  to  absorb  a  single  photon  (6).  Their 
application  of  high-resolution  high-sensitivity  Fourier  transform  infrared  (FTIR)  spectrometry 
has  produced  an  extensive  compilation  of  observations,  covering  most  of  the  600-14000  cm'1 
region,  which  are  unparalleled  in  their  internal  self-consistency.  This  comprehensive  spectral 
survey,  alongside  a  handful  of  classic  prior  solid  parahydrogen  (pH2)  spectroscopic  studies  (7- 
13),  constitute  an  indispensable  resource  to  those  seeking  an  improved  understanding  of 
photodynamics  in  quantum  molecular  solids  —  including  the  growing  number  of  practitioners  of 
Matrix  Isolation  Spectroscopy  (MIS)  in  solid  pH2  matrices. 

The  use  of  pH2  solids  with  very  low  residual  orthohydrogen  (oH2)  concentrations  as  matrix 
hosts  for  chemically  interesting  dopants  originated  in  the  Oka  laboratory  in  Chicago,  and  has 
advanced  in  collaboration  with  Shida  and  Momose  and  coworkers  in  Kyoto  (14-24).  These 
authors  have  reviewed  in  detail  the  unique  properties  of  solid  pH2  as  a  matrix  host,  including  the 
preservation  of  the  spherical  nature  of  the  ground  state  (v  =  0,  J  =  0)  pH2  molecules  in  the  solid 
phase,  and  the  resulting  advantages  to  high  resolution  rovibrational  spectroscopy  (14,16).  Most 
of  the  samples  utilized  in  these  pioneering  studies  were  produced  by  condensation  of  a  room 
temperature  dopant/pH2  gas  mixture  in  an  enclosed  cell  held  near  Vi  the  pH2  triple  point 
temperature  of  Ttp  =  13.8  K.  This  approach  results  in  beautifully  transparent,  centimeters-long, 
polycrystalline  samples  which  are  stable  over  a  relatively  broad  temperature  range. 
Unfortunately,  this  sample  preparation  scheme  also  leads  to  extensive  dopant  aggregation, 
clustering,  and  even  phase  separation. 

A  few  years  ago  we  developed  a  modification  to  the  standard  MIS  substrate-in-vacuum 
sample  deposition  scheme  which  yields  optically  transparent,  millimeters-thick,  doped  pH2 
solids,  and  which  is  compatible  with  most  existing  impurity  doping  methods.  Our  Rapid  Vapor 
Deposition  (RVD)  technique  (25,26),  in  which  volatilized  dopants  are  co-deposited  with  an 
independent  flow  of  precooled  pH2  gas  directly  onto  a  substrate-in-vacuum  held  at  T  «  2  K, 
offers  excellent  dopant  isolation  efficiencies  for  concentrations  up  to  ~  100  ppm  —  ideal  for 
FTIR-based  absorption  diagnostics.  However,  this  improved  isolation  efficiency  comes  at  the 
price  of  a  more  complicated  host  morphology  with  higher  defect  concentrations,  and  a  maximum 


Approved  for  public  release;  distribution  unlimited. 
3 


working  temperature  limited  by  destructive  sample  sublimation  to  T  <  5  K  (23,27).  Despite 
these  shortcomings,  we  are  aware  of  a  number  of  ongoing  efforts  in  independent  laboratories  to 
adapt  the  basic  RVD  scheme  to  a  variety  of  MIS  applications. 

Since  water  vapor  is  the  major  residual  gas  component  in  the  unbaked  o-ring  sealed  high- 
vacuum  systems  typically  employed  in  MIS  studies,  water  is  a  ubiquitous  impurity  in  vapor 
deposited  cryogenic  matrices;  appearing  unintentionally  in  our  RVD  pH2  solids  at  up  to  ~  1  ppm 
levels  (28).  Clearly,  a  priori  knowledge  of  the  spectroscopic  features  attributable  to  the  pH2 
host,  and  to  water  molecules  and  other  likely  contaminants  (29),  is  required  to  avoid  confusion 
while  assigning  spectra  of  new  species  in  solid  pH2,  e.g.  dopant-H2  reaction  products  (30,31). 
Also,  water  clusters  and  water  complexes  with  other  dopants  are  of  fundamental  interest  (28);  the 
process  of  analyzing  their  spectroscopic  signatures  also  begins  with  the  identification  and 
elimination  of  matrix  host  features,  and  more  mundane  water  related  features.  We  hope  that  the 
data  presented  in  this  manuscript  will  facilitate  this  process. 

In  what  follows  we  give  a  brief  review  of  the  rovibrational  spectroscopy  of  matrix  isolated 
JJAJDaO,  and  HDO  molecules,  we  summarize  past,  lessons  learned  aboutthe  microscopic 
structure  of  RVD  pH2  solids,  along  with  a  few  salient  details  of  the  gas  phase  spectroscopy  of 
water  molecules.  We  describe  our  experimental  implementation  of  the  RVD  sample  preparation 
technique,  and  our  FTIR  spectroscopic  diagnostics  covering  the  800-7800  cm'1  range.  We 
present  and  assign  the  direct  IR  absorptions  of  water  monomers  isolated  in  solid  pH2,  including 
data  showing  the  kinetics  of  nuclear  spin  conversion  (NSC)  of  metastable  J  =  1  ortho-H^O 
(0H2O),  para-D20  (pD20),  and  0H2  molecules.  We  discuss  the  observed  matrix  effects  on  the 
water  monomer  spectra,  i.e.  gas-to-matrix  vibrational  band  origin  shifts  and  perturbed  rotational 
bands.  We  assign  a  few  transitions  to  water  molecules  with  a  single  0H2  nearest  neighbor,  i.e. 
“0H2-H2O  complexes,”  and  present  unassigned  spectra  of  larger  (oFfOn-water  clusters  along  with 
a  qualitative  description  of  the  observed  features,  including  NSC  in  oH2-water  complexes  and 
(oH2)n-water  clusters.  Finally,  we  present  and  discuss  IR  features  of  the  pH2  solid  induced  by  the 
presence  of  the  water  dopants,  as  well  as  cooperative  absorption  features  which  arise  when  the 
dopant  and  the  pFb  solid  share  the  energy  of  a  single  IR  photon  (32-34). 
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H.  BACKGROUND 

A.  Matrix  Isolation  Considerations 

As  might  be  expected,  studies  of  trapped  water  monomers  date  back  to  the  origins  of  the 
MIS  field  (35-51).  Rovibrational  spectra  of  water  monomers  in  rare  gas  solids  (RGS)  are  well 
explained  by  assuming  that  the  water  molecules  exist  as  slightly  hindered  rotors  in  single 
substitutional  vacancies  ( 37-41,48-51 ),  except  for  a  few  features  assigned  to  “non-rotating” 
water  molecules  trapped  near  defects  (38,39,50,51).  Water  monomers  trapped  in  molecular 
solids  such  as  N2  and  normal-D2  do  not  rotate  (42-47);  apparently  the  trapping  environments  in 
these  solids  are  significantly  more  anisotropic  than  those  in  RGS.  With  two  important 
exceptions  (i.e.  clustering  with  residual  0H2,  and  the  dopant-induced  and  cooperative  transitions 
associated  with  the  pH2  vibron)  the  spectroscopy  of  water  monomers  in  solid  pH2  is 
understandable  by  comparison  with  this  earlier  RGS  work. 

Since  the  center  of  mass  separation  in  the  gas  phase  H2-H2O  complex  is  R<.  ~  3.6  A  (52),  and 
the  nearest  neighbor  separation  in  solid  pH2  is  Rm  =  3.789  A  (10),  a  water  molecule  should  be 
readily  accommodated  in  a  single  substitutional  vacancy  in  solid  pH2.  Spectra  of  other  dopants 
in  RVD  solid  pH2  have  been  well  explained  by  postulating  that  our  as-deposited  solids  exhibit  a 
mixed  hexagonal-close-packed/face-centered-cubic  (hcp/fcc)  polycrystalline  morphology 
(25,27,29);  the  crystal  field  in  single  substitutional  hep  sites  is  of  D3h  symmetry  (16),  while  fee 
sites  are  of  Oh  symmetry  (27).  However,  the  hep  (...ABCABC...)  and  fee  (...ABABAB...) 
structures  are  only  two  of  the  infinite  number  of  densest  close-packed  structures  derivable  by 
stacking  of  close-packed  planes  (53).  If  only  nearest  neighbor  dopant-pH2  interactions  are 
important,  dopants  in  single  substitutional  sites  in  layers  B  will  exist  in  regions  of  local  hep 
stacking  (...XABCX...,  “hep-like”)  or  local  fee  stacking  (...XABAX...,  “fcc-like”).  Thus,  we 
cannot  rigorously  exclude  structures  based  on  other  stacking  arrangements,  such  as  polytypic  or 
random-stacked  close-packed  (rep)  structures.  Annealing  of  RVD  pH2  samples  at  T  >  4.5  K 
results  in  nearly  complete  and  irreversible  conversion  to  the  thermodynamically  stable  hep 
structure.  We  show  below  that  annealing  of  water  doped  pH2  solids  results  in  only  minor  shifts 
in  peak  positions,  but  often  in  considerable  sharpening  of  the  linewidths. 

B.  Gas  Phase  Water  Molecules 

The  molecular  symmetry  point  group  for  H20  and  D20  molecules  is  C2v,  hence  all  three 
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principal  moments  of  inertia  IA  <  Ib  <  Ic  are  different,  the  three  corresponding  rotational 
constants  A  >  B  >  C,  are  different,  and  these  molecules  exist  as  asymmetric  tops  (54).  The  C2 
symmetry  axis  is  the  intermediate  B  axis,  the  A  axis  lies  in  the  plane  of  the  molecule,  and  the  C 
axis  is  mutually  perpendicular  (55).  FIG.  1  shows  a  few  germane  rovibrational  energy  levels  and 
allowed  IR  transitions  for  gas  phase  H2O  and  D2O  molecules  (37,54).  For  symmetric  tops  the 
rotational  levels  are  labeled  by  the  quantum  numbers  J  and  K,  where  J  is  the  total  rotational 
angular  momentum,  and  K  is  the  constant  component  of  J  about  the  figure  axis  (54).  The  H2O 
and  D2O  asymmetric  rotor  levels  are  labeled  JicaKc,  where  Ka  =  |K|  correlated  to  the  prolate  rotor 
limit  (Ib  =  Ic),  and  Kc  =  |K|  in  the  oblate  rotor  limit  (IA  =  Ib)  (56).  Alternatively,  the  levels  can 
be  labeled  Jx,  where  t  =  Ka  -  Kc,  and  -J  <  x  <  J,  with  x  increasing  with  increasing  level  energy. 
Each  of  these  levels  is  (2J  +  l)-fold  degenerate  in  the  quantum  number  M  arising  from  the 
isotropy  of  free  space. 

The  two  hydrogenic  atoms  are  equivalent  via  the  C2  symmetry  axis,  so  just  as  in  the  case  of 
H2  and  D2  (10)  the  FI20  and  D2O  molecules  each  exist  in  two  distinct  nuclear  spin  (I  =  Ij  +  I2) 
modifications  (54):  para-EEO  (pFEO,  I  =  0),  ortho-FbO  (0H2O,  I  =  1),  orthoriCEO  (oD20, 1  =  0, 
2),  and  para-D20  (pD20, 1  =  1),  where  by  convention  the  modification  with  the  highest 
degeneracy  gi  =  X(2I  +  1)  is  designated  “ortho.”  At  liquid  helium  (IHe)  temperatures  the 
molecules  rapidly  relax  to  the  lowest  rotational  level  available  while  conserving  I:  the  true 
ground  Ooo  state  for  pH^O  and  0D2O,  and  the  metastable  loi  state  for  0H2O  and  pD20. 

The  three  nondegenerate  vibrational  fundamentals  are  labeled  Vi  (symmetric  stretch),  v2 
(bend),  and  v3  (asymmetric  stretch).  Rovibrational  transitions  of  asymmetric  top  molecules  are 
classified  as  Type  A,  Type  B,  or  Type  C,  depending  on  whether  the  transition  dipole  moment  is 
in  the  direction  of  the  least,  intermediate,  or  largest  moment  of  inertia,  respectively  (54).  Thus, 
excitation  of  the  Vi,  v2,  and  v3  fundamentals,  and  all  overtones  and  combination  bands,  are 
further  classified  as  Type  A  bands  if  Av3  is  odd,  or  Type  B  otherwise  (54).  The  operant  IR 
selection  rules  are  summarized  implicitly  in  FIG.  1  as  arrows  showing  the  relevant  dipole 
allowed  transitions.  The  gas  phase  level  energies,  transition  wavenumbers,  and  transition 
probabilities  are  taken  from  the  literature  (57-62). 

The  H  and  D  atoms  in  HDO  are  trivially  nonequivalent,  hence  there  is  no  ortho/para 
distinction,  and  all  HDO  molecules  can  relax  at  low  temperatures  to  the  Ooo  state.  The  molecular 
symmetry  point  group  for  HDO  is  Ci,  and  the  vibrational  modes  change  character  qualitatively, 
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with  the  Vi  mode  resembling  a  localized  O-D  stretch,  and  the  V3  mode  a  localized  O-H  stretch 
(54,63,64)',  so  any  HDO  molecule  can  participate  in  mixed  Type  A  and  Type  B  transitions.  For 
HDO,  the  compilation  of  gas  phase  transition  probabilities  in  the  HITRAN  database  (59)  was 
especially  helpful  for  clarifying  which  transitions  should  be  observable  in  our  experiments. 

III.  EXPERIMENTAL 

Our  experimental  methods  have  been  described  in  detail  before  (25,26).  Briefly:  we 
prepare  our  samples  by  co-deposition  of  independent  flows  of  room  temperature  water  vapor  and 
precooled  pH2  gas  onto  a  BaF2  substrate-in-vacuum  cooled  to  T  *  2  K  in  a  lHe  bath  cryostat.  We 
estimate  the  dopant  concentrations  from  these  gas  flow  rates  (with  a  ±  50  %  uncertainty  in  these 
experiments  due  to  water  adsorption/desorption  within  the  gas  handling  manifold)  and  so  report 
them  here  using  only  one  significant  digit.  In  these  experiments  we  operate  our  ortho/para  H2 
converter  at  T  =  15  K,  yielding  *  100  ppm  residual  0H2  content,  and  at  T  =  19  K  yielding  *  1000 
ppm  residual  0H2.  Our  favorite  200  mmol/hour  pH2  gas  flow  rate  results  in  a  solid  growth  rate 
of*  3  mm/hour.  We  determine  the  sample  thickness  ((55)  within  ±  10  %  from  the  intensities  of 
the  Qi(0)+So(0)  and  Si(0)+S0(0)  double  transitions  (66)  of  the  pFh  host. 

We  measure  IR  absorption  spectra  of  our  samples  across  the  800  to  7800  cm'1  range  using 
an  FTIR  spectrometer  (Bruker  IFS120HR)  equipped  with  glowbar  and  tungsten  lamp  sources, 
Ge:KBr  and  Si:CaF2  beamsplitters,  and  liquid  nitrogen  cooled  HgCdTe  and  InSb  detectors.  To 
accommodate  the  IR  diagnostic,  the  entire  optical  path  is  enclosed  within  a  0.5  m3  polycarbonate 
box  purged  with  a  constant  flow  of  dry  nitrogen  gas.  During  the  FT  process,  we  employ  a 
triangular  apodization  function  which  results  in  a  sin2(x)/x2  instrumental  lineshape  function  (67). 

In  this  study  we  produce  samples  and  collect  IR  data  following  two  different  experimental 
protocols:  (1)  our  standard  approach  of  depositing  a  *  3  mm  thick  sample  over  an  hour  long 
period,  then  collecting  and  co-adding  several  large  interferograms  over  an  additional  hour  to 
yield  high  resolution  (0.005  to  0.01  cm'1)  spectra,  and  (2)  depositing  thinner  samples  in  10  to  15 
minutes  followed  immediately  by  the  rapid  accumulation  of  shorter  interferograms  over  just  2 
minutes  to  yield  low  resolution  (0.05  to  0.1  cm'1)  spectra  favoring  the  observation  of  IR 
transitions  of  metastable  J  =  1  water  molecules.  Individual  sample  preparation  details  are  given 
in  the  figure  captions. 
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rv.  RESULTS  &  DISCUSSION 

A.  H2O  Monomers 

FIG.  2  shows  IR  absorption  spectra  in  the  v2  H2O  region;  a  summary  of  the  observed  peak 
positions,  linewidths,  and  assignments  can  be  found  in  Table  1.  We  make  our  assignments  by 
comparison  with  previous  gas  phase  and  MIS  studies;  the  monomer  peak  positions  in  solid  pH2 
are  uniformly  red  shifted  by  only  «  4  cm'1  from  the  gas  phase  values.  The  0H2O  pH^O  NSC 
process  results  in  diminishing  intensities  for  absorptions  originating  from  the  loi  rotational  level, 
and  growth  of  those  from  the  Ooo  level;  we  will  discuss  the  decay  kinetics  in  greater  detail, 
below.  We  will  also  show  below  that  the  sharp  features  in  FIG.  2  labeled  with  an  “o”  grow  in 
number  and  intensity  with  the  deliberate  addition  of  0H2,  so  we  assign  these  features  to  0H2-H2O 
complexes,  or  possibly  larger  (oFbVHbO  clusters,  identifying  0H2O  and  pH20  where  possible. 
Not  shown  are  spectra  of  the  other  Type  B  bands  of  H20  observed  in  this  study,  2v2  and  vj, 
which  demonstrate  similar  behavior;  these  data  are  also  summarized  in  Table  1. 

Each  of  the  H2O  monomer  features  in  FIG.  2  appears  to  be  a  broad  single  line;  we  observe 
no  lifting  of  the  M  degeneracy  of  the  excited  rotational  levels  by  the  crystal  field.  While  the 
rather  broad  lines  make  it  dangerous  to  dismiss  inhomogeneous  contributions  to  the  linewidths, 
we  note  that  annealing  of  this  sample  results  in  shifting  and  sharpening  of  the  In  <-  Ooo  v2  H20 
line  by  only  ~  0.1  cm'1.  The  nearly  Lorentzian  lineshapes  of  the  1  n  <-  00o  and  212  <-  l0i  v2  H20 
features  are  suggestive  of  the  rovibrational  dephasing  reported  for  CD4  in  solid  pH2  (21). 
Focusing  further  to  the  In  Ooo  feature,  the  measured  annealed  linewidth  of  y  (FWHM)  =  1.6 
cm'1  can  be  related  to  a  dephasing  time  (21)  of  T2=  l/(Ttcy) »  7  ps,  which  is  only  a  few  times  the 
classical  rotational  period  in  the  In  state:  Trot  =  l/c(A  +  C) «  1  ps. 

Annealing  also  leads  to  new  sharp  features  in  the  1595-1600  cm'1  region,  and  to  a  broad 
peak  centered  at  1605  cm'1;  we  tentatively  assign  these  features  to  water  clusters.  Such  water 
cluster  features  will  be  noted  without  significant  discussion  throughout  this  manuscript. 

FIG.  3  shows  IR  absorption  spectra  in  the  v3  H20  region;  again  peak  positions,  widths,  and 
assignments  can  be  found  in  Table  1.  We  discussed  briefly  these  v3  H20  monomer  peaks  in  a 
communication  concerning  water  clusters  in  solid  pH2  (28);  also,  the  l0i  00o  peak  near  3765.5 
cm'1  was  reported  as  the  signature  of  a  H20  photoproduct  in  HI/CCVpHh  samples  (20).  Another 
prior  report  of  H20  peak  positions  in  solid  pH2  lists  line  positions  which  differ  considerably  from 
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ours  (22);  perhaps  due  to  perturbations  by  the  ionic  species  known  to  be  present  in  those 
samples. 

As  in  the  case  of  the  V2  H2O  band,  we  again  observe  single  Lorentzian  peaks  for  each 
monomer  transition,  the  decay  of  0H2O  absorptions  occurs  on  the  ~  1  hour  timescale  of  the 
experiment,  we  note  several  features  attributable  to  0H2-H2O  complexes  or  (0H2VH2O  clusters, 
and  peaks  due  to  water  clusters  appear  upon  annealing  near  3716,  3725,  3738.3,  and  3744.6  cm'1. 
Not  shown  are  spectra  of  the  other  similarly  behaving  Type  A  H2O  bands:  V2+V3,  2V2+V3,  and 
V1+V3;  these  data  are  also  included  in  Table  1 . 

In  comparison  with  the  lower  energy  Vi  and  v2  H2O  bands,  the  two  main  features  of  the  v3 
H2O  band  are  quite  sharp,  suggesting  order-of-magnitude  longer  rovibrational  dephasing  times. 
This  observation  contradicts  the  rough  rule  of  thumb,  proposed  based  on  previously  reported 
spectra  of  dopants  in  solid  pH2,  that  the  lowest  energy  vibrational  mode  should  show  the 
narrowest  lines  for  a  given  dopant  molecule.  This  hypothesis  was  consistent  with  the  picture  of 
intramolecular  vibrational  redistribution,  “IVR”  (68),  playing  an  important  role  in  relaxing  the 
higher  lying  vibrational  modes.  Our  H2O/PH2  observations  indicate  faster  dephasing  of  Type  B 
bands  by  dopant-host  interactions  than  of  Type  A  bands  by  IVR. 

We  note  in  passing  the  momentary  confusion  caused  by  the  near  coincidence  of  the 
Qi(0)+Uo(0)  H2  double  transition  (9),  which  we  observe  at  5320.0  cm'1,  with  the  v2+v3  H20 
combination  band  vibrational  origin. 

B.  D2O  Monomers 

A  summary  of  our  D2O/PH2  data  and  assignments  is  given  in  Table  2.  FIG.  4  shows  the  V2 
D2O  spectral  region  which  is  very  similar  to  the  v2  H2O  spectra,  except:  it  is  more  spectrally 
compact  due  to  the  smaller  rotational  constants  of  D2O,  the  lio  <-  loi  V2  D20  line  is  very  much 
sharper  than  for  V2  H2O,  and  the  decay  of  metastable  J  =  1  D2O  molecules  is  more  than  twice  as 
fast  as  for  1=1  H2O.  FIG.  4  also  shows  another  example  of  an  intruding  solid  hydrogen  line: 
Uo(0)  H2  near  1 167.1  cm'1;  the  increase  in  intensity  of  this  feature  upon  annealing  reflects  the 
irreversible  fee  ->  hep  conversion  of  the  pF^  host.  Clustering  of  D2O  molecules  upon  annealing 
leads  to  new  sharp  peaks  near  1180  and  1 190  cm'1. 

FIG.  5  shows  the  V3  D2O  region  which  again  resembles  a  more  spectrally  compact  version  of 
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the  v3  H20  region,  except  that  the  lines  are  considerably  sharper  in  the  annealed  spectrum.  In 
fact,  upon  annealing  the  10]  <-  000  v3  D20  line  sharpens  to  0.10  cm-1  FWHM,  and  one  of  the 
oH2-oD20  features  actually  split  into  a  triplet  at  2785.05, 2785.15,  and  2785.29  cm'1.  We  also 
note  the  strange  temporal  dependence  of  the  peak  near  2769.3  cm'1,  which  we  assign  to  an  oH2- 
pD20  complex,  but  which  does  not  follow  the  decay  kinetics  of  pD20.  This  feature  is  only 
present  in  D20  doped  samples,  so  it  is  not  a  transition  of  the  pH2  host;  we  speculate  that  it  may 
overlap  a  (D20)n  cluster  line,  such  as  the  ones  apparent  near  2751, 2764,  and  2776  cm'1. 

C.  Nuclear  Spin  Conversion  of  Isolated  H2O  &  D2O  Monomers 

FIG.  6  plots  the  integrated  intensities  of  a  few  transitions  of  J  =  1  H20  and  D20  monomers 
vs.  time  after  sample  deposition.  We  present  only  those  peaks  which  do  not  overlap  water  cluster 
features  and  so  decay  cleanly  to  zero  baseline.  At  T  =  2.4  K,  the  oH20  population  decays  as  a 
single  exponential  with  a  lifetime  of  1900  ±  100  s,  the  pD20  population  decays  as  a  single 
exponential  with  a  lifetime  of  860  ±  50  s.  Not  shown  are  plots  of  the  growth  of  the  J  =  0  water 
monomer  transitions  which  show  the  same  single  exponential  first  order  kinetics.  We  did  not 
investigate  the  temperature  dependence  of  these  NSC  rates,  nor  any  possible  photo-effects  due  to 
room  lights  or  the  IR  radiation  from  the  spectrometer. 

Previous  studies  of  NSC  of  H20  in  solid  Ar  matrices  {39,41)  reported  much  longer 
conversion  times  than  ours,  unless  an  02  dopant  was  deliberately  added  at  0.1  - 1  % 
concentrations.  Those  authors  concluded  that  interactions  with  the  paramagnetic  02  molecules 
accelerate  the  rate  of  decay  of  oH20;  we  do  not  believe  that  02  impurities  play  any  role  in  our 
experiments.  Based  on  the  measured  composition  of  residual  gases  in  our  cryostat  prior  to 
cooldown,  we  expect  02  contamination  to  be  smaller  than  H20  contamination  by  at  least  a  factor 
of  10;  putting  an  upper  bound  on  trapped  02  concentrations  of-  100  ppb.  Furthermore,  we  do 
not  see  any  spectroscopic  evidence  for  02-H20  complexes,  while  we  clearly  do  so  for  oH2-H20 
complexes  {vide  infra). 

Our  results  more  closely  resemble  those  reported  for  NSC  of  CH4  molecules  in  solid  pH2 
where  clean  first  order  kinetics  were  also  observed  {19).  This  study  also  documented  a  curious 
dependence  of  the  NSC  rate  constant  on  the  concentration  of  oH2  molecules:  essentially  no 
dependence  for  oH2  concentrations  between  0.05  and  2  %,  then  a  linearly  increasing  rate 
constant  with  oH2  concentrations  of  10  to  75  %.  The  authors  conclude  that  the  increased  NSC 
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rate  is  due  to  magnetic  dipole-dipole  interactions  between  nearest  neighbor  0H2  and  methane 
molecules,  but  the  mechanism  for  NSC  at  the  lower  0H2  concentrations  remains  unclear  (19). 
More  recent  unpublished  work  by  this  group  has  demonstrated  the  strong  temperature 
dependence  of  NSC  rates  in  low  0H2  concentration  samples,  indicating  that  interactions  between 
the  J  =  1  water  dopants  and  the  non-magnetic  pH2  host  are  capable  of  promoting  the  NSC 
process  (69). 

Since  the  0H2  concentrations  in  the  samples  depicted  in  FIGS.  2-6  are  only  100  ppm,  and 
with  the  benefit  of  clear  spectroscopic  signatures  of  isolated  water  monomers  as  distinct  from 
oH2-water  complexes,  we  assert  that  our  reported  NSC  lifetimes  are  indeed  for  0H2O  and  PD2O 
molecules  having  no  0H2  nearest  neighbors.  These  data  may  provide  a  valuable  test  for  models 
of  NSC  invoking  “non-magnetic”  mechanisms  (70-72).  We  will  discuss  the  NSC  of  (oFkV 
0H2O  clusters  further  below. 

D.  HDO  Monomers 

A  summary  of  HDO/pFt  peak  positions,  linewidths,  and  assignments  appears  in  Table  3. 
FIG.  7  shows  the  v2  HDO  region  of  the  IR  spectrum  of  a  mixed  H2O/HDO/D2O/PH2  sample 
which  is  dominated  by  two  monomer  peaks,  which  we  assign  to  the  l0i  4-  Ooo  and  In  4-  O0o 
transitions  based  on  their  2  to  3  cm'1  gas-to-matrix  red  shifts;  these  shifts  being  intermediate 
between  those  observed  for  v2  H2O  and  V2  D2O  features.  Since  HDO  does  not  exist  in  ortho/para 
modifications,  at  T  =  2.4  K  we  expect  99.99  %  of  the  HDO  molecules  to  rapidly  relax  to  the 
ground  Ooo  (v  =  0)  level.  Using  the  gas  phase  rotational  excitation  energy  for  the  l0i  (v  =  0)  level 
of «  15.5  cm'1,  we  predict  a  100  x  increase  in  thermal  population  for  warming  from  T  =  2.4  K  to 
T  =  4.8  K.  Indeed,  warming  this  sample  to  T  =  4.8  K  results  in  the  appearance  of  the  Ooo  <—  loi 
peak,  again  at  a  «  3  cm'1  gas-to-matrix  red  shift.  Despite  comparable  or  even  larger  transition 
strengths  (59),  we  observe  no  other  V2  band  absorptions  from  the  10]  level;  these  bands  are  either 
too  broad  or  overlapped  by  other  absorption  features.  Repeated  temperature  cycling  shows  these 
spectral  changes  to  be  completely  reversible,  with  the  exceptions  of  the  dramatic  sharpening  of 
the  absorption  lines  and  the  appearance  of  water  cluster  features  near  1399  and  1407  cm'1.  The 
0H2-FEDO  complex  peak  at  1413.28  cm'1  resolves  into  a  doublet  at  1413.20  and  1413.29  cm'1 
upon  annealing. 

FIG.  8.  shows  the  V3  HDO  region  which  is  initially  overlapped  by  V3  0H2O  and  0H2-0H2O 
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complex  absorptions.  We  again  observe  a  transition  from  the  thermally  populated  loi  level  upon 
warming  the  sample  to  T  =  4.8  K.  Annealing  again  results  in  dramatic  sharpening  of  the 
absorption  lines,  with  the  oH2-HDO  multiplet  near  3705.5  cm'1  resolving  into  a  triplet  of  sharp 
(0.03  cm'1  FWHM)  peaks  at  3705.39,  3705.44,  and  3705.53  cm'1,  and  new  water  cluster  features 
appearing  near  3700  cm’1. 

The  annealed  linewidths  for  HDO  monomers  are  the  sharpest  we  have  encountered  for  any 
of  the  water  monomers,  perhaps  due  to  the  qualitative  differences  in  vibrational  mode  character 
and  IVR  pathways  for  the  mixed  isotopomer  (54,63,64)-  FIG.  9  shows  the  effects  of  annealing 
on  the  loi  4-  00o  Vi  HDO  line,  which  splits  into  two  peaks:  a  main  Lorentzian  line  at  2729.995 
cm'1,  and  a  small  satellite  at  2730.085  cm'1;  the  0.03  cm'1  linewidth  of  the  main  peak  yields  a 
dephasing  time  of  T2  «  0.3  ns.  We  suggest  that  the  main  peak  is  due  to  HDO  monomers  in  hcp- 
like  trapping  sites,  and  the  satellite  to  HDO  monomers  in  residual  fcc-like  trapping  sites.  The 
0.09  cm'1  splitting  between  these  peaks  is  consistent  with  our  previous  estimate  of «  0.1  cm'1 
residual  inhomogeneous  linebroadening  of  the  other  water  monomer  transitions. 

E.  Matrix  Effects  on  Monomer  Spectra 
1.  Gas-Matrix  Shifts 

FIG.  10  shows  a  correlation  plot  of  water  monomer  gas-to-matrix  red  shifts  vs.  rdvibrational 
excitation,  summarizing  these  data  from  Tables  1,  2,  and  3.  The  correlation  is  monotonic  and 
roughly  linear,  however  we  can  readily  discern  two  distinct  subtrends:  the  slope  for  v2  excitation 
is  (Egas-Ematrix)/Egas «  0.002,  whereas  the  slope  for  Vi  or  V3  excitation  is  »  0.0035;  these  trends  are 
maintained  for  overtones  and  combination  bands.  We  do  not  understand  these  trends,  and 
present  our  data  and  reasoning  here  to  highlight  this  unresolved  question. 

We  initially  expected  these  trends  to  be  (at  least  qualitatively)  easily  explicable  using  the 
standard  treatment  of  matrix  shifts  of  the  vibrational  band  origins,  which  assumes  the  separation 
of  dopant/host  intermolecular  interactions  into  a  sum  of  electrostatic  +  induction  +  dispersion  + 
repulsive  terms  (73): 

V  —  Velec  Vinc)  +  Vdisp  +  Vrep  [1]. 

The  gas-to-matrix  band  origin  shifts  are  then  calculated  perturbatively  by  considering  the 
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changes  to  each  of  these  terms  upon  vibrational  excitation  of  the  dopant  molecule.  Since  1  =  0 
pH2  molecules  are  spherically  symmetric  and  have  no  rotationally  averaged  permanent 
electrostatic  multipole  moments,  we  can  safely  neglect  the  electrostatic  interaction  term.  The 
repulsive  interactions  are  the  most  difficult  to  evaluate,  and  are  likely  to  increase  in  magnitude 
upon  vibrational  excitation  of  the  dopant,  contributing  to  spectral  blue  shifts;  we  will  tentatively 
ignore  them  for  the  moment. 

The  electrostatic  induction  solvation  contribution  to  the  vibrational  band  origin  matrix  shift 
is  predicted  to  be  proportional  to  the  difference  in  the  squares  of  the  dopant’s  molecular  electric 
dipole  moment  upon  vibrational  excitation  (75): 

Ay  Vind  oc  Av  (p  )  [2] 

where:  Av  (p2)  =  [p2]v  -  [p2]o,  and  p  is  the  rotationless  molecular  electric  dipole  moment.  The 
applicability  of  this  expression  to  individual  rovibrational  transitions  is  suspect,  due  to  rotational 
averaging  of  the  dopant  molecule’s  dipole  moment.  However,  our  observation  of  only  slightly 
matrix  perturbed  rotational  band  structures  (vide  infra )  means  that  the  individual  rovibrational 
transitions  indeed  shift  with  the  vibrational  band  centers.  The  known  vibrational  dependence  of 
p  (74, 75)  predicts  red  shifts  for  vi  and  V3  excitation:  Av  vind  (vi)  «  Av  vmd  (V3)  and  a  blue  shift 
(i.e.  a  decrease  in  |p|)  upon  v2  excitation,  so  induction  interactions  alone  cannot  explain  the 
observed  trends  in  matrix  shifts.  Furthermore,  this  discrepancy  cannot  be  explained  by  rotational 
averaging  effects,  as  the  vi  and  v2  bands  share  the  same  rotational  structure. 

The  vibrational  band  matrix  shift  due  to  changes  in  dispersion  energies  is  proportional  to  the 
change  in  the  dopant’s  molecular  polarizability  upon  excitation  (75): 

VdiSp  OC  Ay  (cc)  [3] 

where:  Av  (a)  =  [a]v  -  [a]o,  and  a  is  the  mean  H20  polarizability.  The  vibrational  dependence 
of  the  mean  polarizability  of  H20  (76)  predicts  Av  vdisp  (vi) «  Av  vdisp  (v2)  <  Av  vdisp  (v3),  which 
again  incorrectly  places  the  contribution  from  excitation  of  the  Vj  symmetric  stretch  mode  well 
below  that  for  the  V3  asymmetric  stretch  mode. 
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We  see  that  neither  the  induction  or  dispersion  contributions,  either  singly  or  in 
combination,  can  explain  the  details  of  the  observed  matrix  shifts.  A  proper  calculation  will 
likely  involve  a  combination  of  all  the  interactions  in  the  standard  model,  including  the  repulsive 
interactions  and  perhaps  even  more  subtle  effects  such  as  hydrogen-bonding  H2-H2O 
interactions,  which  is  beyond  the  scope  of  the  present  study. 

2.  Rotational  Analysis 

Since  we  only  observe  IR  transitions  originating  from  the  Ooo  and  loi  rotational  levels  in  the 
vibrational  ground  state  for  H20,  D20,  and  HDO  monomers,  we  have  very  limited  information 
with  which  to  conduct  a  rotational  analysis.  We  did  observe  a  handful  of  transitions  originating 
in  a  common  loi  lower  state,  which  yields  a  few  rotational  level  splittings  in  the  vibrationally 
excited  states;  these  data  are  summarized  in  Table  4.  We  employ  the  Mecke  sum  rules  for  the 
rigid  asymmetric  top  ( 54)  to  recast  the  observed  energy  splittings  as  sums  of  rotational  constants. 
We  note  a  small  (and  given  the  large  widths  of  the  observed  absorption  features  perhaps  not 
significant)  gas-to-matrix  increase  in  the  splittings  between  the  2 12  and  ho  rotational  levels  (i.e. 
in  the  C  rotational  constant)  in  the  excited  v2  and  V]  vibrational  states  for  H20;  all  other 
measured  splittings  and  rotational  constant  sums  showed  modest  2  to  5  %  decreases. 

We  were  also  able  to  bracket  the  vibrational  origin  for  a  number  of  transitions.  We  define 
an  energy  difference  equal  to  the  sum  of  the  rotational  energies  of  the  loi  level  in  the  ground  and 
vibrationally  excited  states: 

AEasym  =  v  [hi  Ooo]  ■  V  [Ooo  hi]  [4] 

where  v  signifies  the  transition  wavenumber  of  the  Av  =  +1  vibrational  excitation  accompanied 
by  the  bracketed  rotational  transition.  We  can  adapt  the  Mecke  sum  rules  to  estimate 
vibrationally  averaged  rotational  constants: 

AEasym  -  (B  +  C)v=l  +  (B  +  C)v=o  =  2  (B  +  C)avg  [5]. 

These  results  are  summarized  in  Table  5  which  shows  a  consistent  w  3  %  gas-to-matrix 
reduction,  again  validating  our  picture  of  water  monomers  in  solid  pH2  existing  as  very  slightly 
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hindered  rotors. 


F.  0H2-0H2O  Complexes  and  (ofyjn-water  Clusters 

We  confirm  the  identification  of  the  features  labeled  with  an  “o”  in  FIGS.  2-8  to  absorptions 
of  0H2-0H2O  complexes  and  (oH2)n-H20  clusters  by  deliberately  doping  some  samples  with  « 
1000  ppm  0H2,  which  we  accomplish  by  operating  the  ortho/para  converter  at  T  =  19  K.  Spectra 
of  the  V2  H2O  region  are  shown  in  FIG.  11,  and  spectra  of  the  V3  H2O  region  are  shown  in  FIG. 
12.  The  features  appearing  at  the  lowest  0H2  concentrations  are  tentatively  attributed  in  Tables 
1-3  to  0H2-H2O  complexes,  although  some  of  the  weaker  features  may  well  be  due  to  water 
clustering  with  multiple  0H2  molecules. 

We  note  in  passing  that  (assuming  preservation  of  natural  isotopic  abundances)  these  pH2 
solids  also  contain  »  300  ppm  of  a  benign  and  immobile  HD  impurity,  which  rapidly  relaxes  to  J 
=  0  upon  cooling,  making  it  virtually  indistinguishable  from  the  pH2  host  molecules.  In  contrast, 
0H2  molecules  have  a  nonzero  rotationally  averaged  electrostatic  quadrupole  moment  (10),  so 
they  can  interact  more  strongly  with  dopant  molecules  than  can  the  spherically  symmetric  J  =  0 
pH2  molecules,  resulting  in  some  clustering  during  the  deposition  process.  Also,  0H2  molecules 
are  mobile  through  the  pH2  solid  by  a  process  known  as  “quantum  rotational  diffusion”  in  which 
an  oH2-pH2  nearest  neighbor  pair  converts  into  a  pH2-oH2  pair  (8,10). 

The  cleanest  examples  of  the  spectroscopic  consequences  of  oH2-dopant  clustering  are 
found  in  a  recent  report  on  (oH2)n-CH3F  clusters  in  solid  pH2  (77).  CH3F  does  not  rotate  in  solid 
pH2,  so  the  rovibrational  spectrum  of  the  v3  C-F  stretching  mode  collapses  to  a  single  line  at  the 
vibrational  band  origin.  The  quadrupolar  oH2  molecules  bind  strongly  to  the  electrostatic  dipole 
of  the  non-rotating  CH3F  molecule,  and  each  of  the  1 3  observed  spectral  features  can  be 
associated  with  a  CH3F  molecule  with  from  0  to  12  oH2  nearest  neighbors. 

At  first  glance  the  spectra  of  (oH2)n-H20  clusters  in  the  v2  H20  region  appear  to  tell  a 
similar  story.  The  inset  to  FIG.  1 1  shows  the  detail  near  the  v2  H20  vibrational  origin  where  the 
intensity  and  number  of  features  increase  with  addition  of  oH2,  however  not  in  a  neat  countable 
fashion  as  for  (oH2)n-CH3F.  The  sharp  features  near  the  vibrational  origin  are  indeed  evocative 
of  the  “non-rotating”  water  absorptions  reported  for  water  molecules  trapped  near  defects  in 
RGS  (42-47).  This  “strictly  non-rotating”  interpretation  requires  that  the  anisotropy  of  the 
trapping  environment  for  an  H20  molecule  adjacent  to  even  a  single  oH2  molecule  be  large 
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enough  to  completely  quench  the  H2O  rotation,  leaving  it  to  oscillate  about  some  space-fixed 
equilibrium  orientation  (75).  In  this  limit,  angular  momentum  is  conserved  at  the  level  of  the 
solid  as  a  whole,  and  the  dopant  angular  momentum  quantum  numbers  become  moot.  Indeed, 
the  total  integrated  intensity  for  the  0H2-H2O  complex  peaks  in  the  1593  to  1595  cm'1  region  in 
FIG.  1 1(a)  decays  much  more  slowly  (only  about  10  %  in  1  hour)  than  do  the  0H2O  monomer 
features,  and  obviously  do  not  grow  in  the  manner  of  PH2O  monomers,  suggesting  either  (1)  that 
the  0H2-0H2O  and  0H2-PH2O  spectroscopic  features  overlap,  or  perhaps  (2)  that  there  is  no 
surviving  distinction  between  0H2-0H2O  and  0H2-PH2O  complexes. 

This  appealingly  simple  interpretation  is  refuted  by  the  spectra  in  FIG.  12,  which  show  the 
qualitatively  different  behavior  of  the  V3  H2O  band  upon  addition  of  0H2.  In  this  case  we  see 
very  little  IR  activity  at  the  vibrational  band  origin,  and  instead  observe  progressions  of  peaks 
building  inwards  from  each  of  the  H2O  monomer  loi  <—  Ooo  and  Ooo  <—  loi  transitions.  In  fact,  the 
features  at  3725.46  and  3758.46  cm'1  appear  at  low  0H2  concentrations,  and  display  the  same 
NSC  kinetics  as  the  nearby  water  monomer  transitions  (vide  infra),  suggesting  that  they 
correspond  to  absorptions  of  0H2-0H2O  and  0H2-PH2O  complexes,  respectively. 

We  can  harmonize  the  observations  for  the  V2  H2O  and  V3  H2O  bands  by  incorporating 
elements  of  what  is  known  about  gas  phase  H2-H2O  complexes,  taking  into  account  a  couple  of 
obvious  and  subtle  condensed  phase  effects.  The  ground  state  of  the  gas  phase  PH2-H2O 
complex  resembles  the  Ar-F^O  complex  and  is  best  described  as  a  pseudo-diatomic  with  very 
slightly  hindered  rotation  of  the  H2O  subunit  (52).  The  ground  state  of  the  gas  phase  0H2-H2O 
complex  exhibits  a  nontrivial  barrier  to  in-plane  rotation  of  the  H2O  subunit,  with  an  equilibrium 
geometry  resembling  the  minimum  energy  “H2O  proton  acceptor”  C2V  configuration  in  which  the 
H2  bond  points  along  the  C2  symmetry  axis  towards  the  O  atom  (52).  However,  large  amplitude 
internal  rotor  dynamics  of  the  H2  and  H2O  subunits  results  in:  (1)  a  ground  state  0H2-PH2O 
wavefunction  with  the  H2  molecule  being  only  ~  8  times  more  likely  to  be  on  the  O  side  than  the 
H  side  of  H2O,  and  (2)  a  node  at  the  C2V  geometry  for  the  lowest  energy  0H2-0H2O  state, 
precluding  an  accurate  assignment  of  the  gas  phase  spectra  using  a  semirigid  C2V  asymmetric  top 
model  (52). 

The  most  obvious  consequence  of  the  matrix  “cage  effect”  on  an  0H2-H2O  complex  trapped 
in  solid  pH2  is  quenching  of  the  end-over-end  rotation  of  the  complex  relative  to  the  crystal 
lattice;  although  “internal”  rotation  and  reorientation  of  the  0H2  and  H2O  subunits  may  not  be 
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strongly  perturbed.  Additionally,  axial  confinement  of  the  0H2  and  H2O  subunits  could  reduce 
their  average  separation,  leading  to  increased  anisotropic  interactions. 

A  more  subtle  matrix  effect  is  the  partial  cancellation  of  0H2-H2O  anisotropic  interactions  by 
PH2-H2O  interactions  in  the  highly  symmetrical  trapping  environment  (15).  This  is  imperfectly 
illustrated  by  considering  the  net  anisotropic  potential  for  rotation  of  H2O  within  a  linear  0H2- 
H2O-PH2  construct.  More  formally,  the  anisotropic  potential  experienced  by  an  H20  molecule  in 
an  0H2-H2O  complex  trapped  in  solid  pH2  can  be  written  in  analogy  to  Ref.  (15)  as: 

Vaniso  =  VoH2(Rl,Ql,®)  +  Xi>i  VpH2(Ri>^)  [6] 

where  the  H2O  molecule  occupies  the  coordinate  system  origin,  the  0H2  molecule  sits  at  position 
Ri  with  orientation  Qi  in  the  crystal  frame,  ©  is  the  orientation  of  the  H2O  molecule  relative  to 
the  0H2-H2O  pair  axis  and  Q  the  H2O  orientation  relative  to  the  crystal  frame,  and  the  remaining 
pH2  neighbors  sit  at  positions  R;.  The  partial  cancellation  effect  is  most  clearly  represented  by 
adding  and  subtracting  the  interactions  of  a  hypothetical  extra  pfk  molecule  at  the  oH2  site  Ri 
(15): 

Vaniso  =  VoffiCRlA,©)  '  VpH2(Rl,0j)  +  ^ .  VpH2  (R,,T1  )  [7] 

or: 

Vaniso  —  A V oH2-pH2(Rl  1 ,®)  "b  VcryStal(f2)  [8] 

where  AV0h2-ph2  is  the  difference  in  anisotropic  interactions  between  an  0H2-H2O  and  a  PH2-H2O 
pair,  and  Vcrystai  is  the  pure  crystal  field  term. 

Because  of  the  cage  effect,  the  partial  cancellation  effect,  and  the  pure  crystal  field  term,  we 
cannot  adopt  unmodified  the  rotational  dynamics  and  spectroscopic  features  of  the  gas  phase 
0H2-H2O  complexes  to  explain  our  condensed  phase  situation.  In  the  language  of  Ref.  (52), 
these  additional  interactions  would  modify  the  magnitude  of  the  angular  anisotropy  parameter,  Q, 
connecting  the  free  H2O  rotor  (C,  =  0)  and  strongly  hindered  rotor  (C,  =  1)  limits.  It  is  difficult  to 
judge  whether  the  net  effect  would  be  to  increase  or  decrease  C,;  the  required  condensed  phase 
simulations  are  beyond  the  scope  of  this  study,  so  we  abandon  any  quantitative  assignment  of  the 
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oH^-water  complex  spectra. 

However,  by  assuming  a  semirigid  C2V  structure  for  the  ground  state  0H2-H2O  complex  in 
the  solid  (i.e.  C,  ~  1),  we  can  qualitatively  understand  the  0H2-H2O  complex  features  in  the  V2  and 
v3  regions  as  parallel  and  perpendicular  bands,  respectively,  of  an  asymmetric  top  near  the 
prolate  symmetric  top  limit  (54).  In  this  description  the  C2  symmetry  axis  now  coincides  with 
the  least  moment  of  inertia  and  is  relabeled  as  the  A  axis,  while  rotations  about  the  new  B  and  C 
axes  are  rendered  infeasible  by  the  matrix  surroundings.  Thus,  the  rotational  constant  for  the 
complex  around  the  C2  symmetry  axis  should  be  comparable  to  the  B  constant  of  the  water 
monomer,  A0H2-H20  «  Bh20,  and  the  new  B  and  C  rotational  constants  for  the  complex  are 
effectively  zero.  The  rotational  term  energies  for  the  symmetric  top  are  given  by  (54): 

E(J,K)  =  BJ(J  + 1)  +  (A  -  B)K2  AK2  asB->0  [9], 

Thus,  we  see  that  the  features  bunching  near  the  V2  H20  band  origin  arise  from  the  AK  =  0 
selection  rule  for  parallel  bands  which  causes  the  0H2-PH2O  and  oH2-oH20  peaks  to  overlap, 
while  the  more  widely  separated  oH2-pH20  and  oH2-oH20  features  in  the  v3  H2O  region  arise 
from  the  AK  =  ±  1  selection  rule  for  perpendicular  bands. 

With  this  insight  we  can  revisit  the  analysis  of  “vibrationally  averaged  rotational  constants” 
presented  above  in  Eqs.  [4]  and  [5]  for  the  water  monomers,  and  apply  it  to  the  main  oH2-water 
complex  peaks  straddling  the  v3  H2O  and  v3  D2O  band  origins.  Using  symmetric  top  quantum 
numbers  (J,K)  we  assign  the  oH2-water  complex  peaks  just  to  the  blue  of  the  00o  4-  loi  water 
monomer  transitions  as  (0,0)  4-  (1,1),  and  those  just  to  the  red  of  the  l0i  4-  00o  monomer 
transition  as  (1,1)  4-  (0,0).  We  note  that  the  relabeling  of  the  rotational  axes  spoils  the  expected 
correlation  between  the  free  rotor  Ka  and  the  symmetric  top  K  quantum  numbers.  We  define 
another  energy  difference: 

AEsym  ■  v  [(1,1)  4-  (0,0)]  -  v  [(0,0)  4-  (1,1)]  [10] 

where  v  again  signifies  the  transition  wavenumber.  From  Eq.  [9]  we  can  relate  this  energy 
difference  to  the  A  rotational  constants  of  the  oH2-water  complex  in  the  ground  and  excited 
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vibrational  states: 


AESym  —  (A)v=l  +  (A)v=0  2(A)avg  [11]- 

From  the  data  in  Tables  1  and  2  we  find  Aavg(v3  oH2-H20)  =  16.5  cm'1  as  compared  with  the 
monomer  Bavg(H20)  «  14.5  cm'1,  and  Aavg(v3  oH2-D20)  =  8.0  cm'1  as  compared  with  the 
monomer  Bavg(D20)  «  7.3  cm'1.  The  reduction  of  the  IA  moment  of  inertia  in  the  complex 
relative  to  Ib  for  the  water  monomers  is  a  consequence  of  the  semi-rigid  nature  of  the  complex; 
the  H(D)  atoms  are  not  always  fully  extended  away  from  the  C2  symmetry  axis.  We  judge  this 
level  of  (dis)agreement  as  acceptable,  given  the  crude  nature  of  our  model. 

The  small  (~  0.1  cm'1)  splittings  in  these  oH2-water  complex  features  observed  in  annealed 
samples  are  probably  due  to  a  combination  of  the  oH2  molecules  occupying  in-plane  (ip)  and 
out-of-plane  (oop)  positions  relative  to  the  H20  subunit  in  the  hep  pH2  lattice  (75),  and  to  lifting 
of  the  M-degeneracy  of  the  symmetric  top  rotational  levels  by  the  crystal  field.  A  more  detailed 
assignment  of  these  splittings,  and  of  the  minor  features  attributed  to  oH2-water  complexes  in  the 
Tables,  will  require  quantitative  modeling  of  the  cage  effect,  and  the  partial  cancellation  and  pure 
crystal  field  terms  in  Eq.  [7]. 

We  can  qualitatively  explain  the  spectral  shifts  of  the  other  (oH2)n-H20  cluster  features  in 
the  v3  region  by  a  simple  electrostatic  argument.  The  strength  of  the  binding  between  oH2  and 
H20  molecules  depends  on  the  degree  of  mutual  orientation  between  the  oH2  and  H20  subunits, 
which  increases  with  the  H20  subunit  angular  momentum.  Thus,  the  features  in  the  3725  to 
3735  cm'1  region  are  due  to  (oH2)n-oH20  clusters,  and  these  absorptions  are  blue  shifted  due  to 
electrostatic  stabilization  of  the  J  =  1,  v  =  0  initial  level.  Conversely,  those  in  the  3745  to  3760 
cm'1  region  are  due  to  (oH2)n-pH20  clusters,  and  these  absorptions  are  red  shifted  due  to 
stabilization  of  the  J  =  1,  v  =  1  final  level. 

G.  Nuclear  Spin  Conversion  in  oH.2-water  Complexes  and  (oH2)„-water  Clusters 

Surprisingly,  NSC  of  metastable  oH20  or  pD20  molecules  is  not  accelerated  appreciably  by 
complexation  with  a  single  oH2  molecule.  This  argues  against  the  importance  of  pairwise 
magnetic  dipole-dipole  interactions  in  the  water  NSC  process,  which  illuminates  the  observed 
independence  of  the  rate  of  NSC  for  CH4  in  solid  pH2  for  oH2  concentrations  below  2  %  (19).  A 
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correlation  plot  of  the  integrated  intensities  of  the  Ooo  <—  loi  V3  H2O  monomer  peak  at  3719.79 
cm'1  and  of  the  (0,0)  <—  (1,1)  0H2-0H2O  complex  peak  at  3725.46  cm'1  shows  a  strong  linear 
correlation  (r2  =  0.994)  over  the  two  hour  duration  of  the  experiment  depicted  in  FIG.  6. 
Furthermore,  the  decay  of  the  (0,0)  <—  (1,1)  0H2-0H2O  peak  at  3725.46  cm'1  is  matched  by  the 
growth  of  the  (1,1)  <-  (0,0)  0H2-PH2O  peak  at  3758.46  cm'1;  the  sum  of  the  integrated  intensities 
of  these  two  features  decays  slowly  by  about  10  %  per  hour,  matching  the  slow  decay  of  the 
overlapped  0H2-0H2O  and  0H2-PH2O  peaks  in  the  1593  to  1595  cm'1  section  of  the  V2  H2O 
region.  This  suggests  that  0H2-0H2O  complexes  convert  directly  into  0H2-PH2O  complexes,  and 
that  absorption  features  for  both  species  have  nearly  the  same  transition  strength.  We  suppose 
that  the  10  %  per  hour  decrease  in  total  0H2-H2O  intensity  is  due  to  NSC  of  0H2  complexed  with 
H20. 

Even  more  surprising  is  the  observation  that  the  integrated  intensity  of  the  bunch  of  (oFEV 
0H2O  cluster  features  in  the  3725  to  3735  cm'1  region  shown  in  FIG.  12(b)  decays  more  slowly 
than  do  the  0H2O  monomer  absorptions!  This  is  just  the  opposite  of  what  we  expected  after 
learning  of  the  increased  rate  of  NSC  for  CH4  in  solid  pH2  for  oH2  concentrations  above  10  % 
(19).  Unfortunately,  we  did  not  investigate  the  kinetics  of  this  process  in  detail;  however,  the 
integrated  intensity  of  the  (oH2)n-oH20  cluster  features  in  FIG.  12(b)  decreased  by  only  50  % 
after  the  first  hour,  and  by  only  70  %  after  five  hours.  We  speculate  that  somehow  the  NSC 
processes  for  the  0H2  and  0H2O  species  are  linked,  with  the  longer  lived  oH2  species  providing  a 
“spin  reservoir”  to  replenish  the  oH20  population. 

H.  Dopant-Induced  IR  Activity  and  Cooperative  Transitions 

Dopant  induced  IR  activity  and  dopant-host  cooperative  absorptions  in  RVD  pFE  solids 
containing  neutral  chemically  interesting  dopants  have  been  reported  and  discussed  previously 
(25,32-34, 77).  FIG.  13.  shows  the  dependence  of  H20-induced  IR  activity  on  NSC  of  the  H20 
dopant;  the  peak  positions  and  assignments  are  summarized  in  Table  6.  The  sharp  Qi(0)  H2 
transition  at  4149.03  cm'1  induced  by  pfUO  is  characteristic  of  a  spherically  symmetric  J  =  0 
perturber,  in  which  case  the  only  possible  induction  mechanism  is  short  range  overlap 
interactions  (32-34).  The  sharpness  of  this  line  indicates  that  dopant-host  interactions  have 
resulted!  in  strong  (several  cm'1)  shifting  of  the  vibrations  of  the  nearby  participating  pH2 
molecules,  detuning  them  from  the  delocalized  vibron  band  of  the  solid  (34).  The  broader,  red- 
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degraded  Qi(0)  H2  feature  peaking  near  4149.8  cm'1  induced  by  0H2O  is  characteristic  of  a 
perturber  possessing  nonzero  rotationally  averaged  electrostatic  multipole  moments,  enabling 
long  range  dopant-pH2  interactions.  The  asymmetric  “vibron-hopping”  lineshape  is  reminiscent 
of  the  familiar  0H2  electric-quadrupole-induced  Qi(0)  H2  line  near  4153.1  cm'1  (11,32-34). 

Working  from  these  origins,  we  see  a  number  of  cooperative  transitions  in  which  a  pure  H2 
vibrational  transition  is  accompanied  by  a  pure  rotational  transition  of  the  dopant  molecule.  This 
has  the  net  effect  of  “upshifting”  the  dopant’s  rotational  spectrum  from  the  far-IR  to  the  mid-IR, 
although  the  rotational  features  observed  in  this  manner  are  for  a  trapped  water  molecule 
additionally  perturbed  by  the  pH2  vibrational  excitation.  The  lio  <-  loi  0H2O  rotational 
transition  is  the  most  strongly  perturbed,  showing  an  18  %  decrease  from  the  gas  phase;  the  same 
transition  also  shows  the  broadest  linewidth  in  the  V2  H2O  rovibrational  band.  Otherwise  we  see 
2  to  5  %  reductions  of  gas  phase  rotational  spacings,  and  2  to  3  cm'1  linewidths,  all  comparable 
to  the  values  reported  above  in  Tables  4  and  5  for  the  direct  IR  rovibrational  Type  B  transitions; 
thus,  rotational  dephasing  (as  distinct  from  rovibrational  dephasing)  is  the  dominant 
linebroadening  mechanism  for  all  these  transitions.  All  we  can  say  at  present  about  the  weak 
features  near  4140  cm'1  is  that  they  are  due  to  the  presence  of  0H2-H2O  complexes. 

FIG.  14.  shows  the  dependence  of  D20-induced  IR  activity  on  NSC  of  the  D2O  dopant,  and 
also  includes  some  HDO  related  features  which  remain  constant  for  the  duration  of  the 
experiment.  All  of  these  observations,  and  their  consequences,  are  very  similar  to  those  for  the 
H2O  dopant,  allowing  for  the  smaller  rotational  constants  of  HDO  and  D2O.  The  HDO  induced 
Qi(0)  H2  peak  is  probably  overlapped  by  the  corresponding  D2O  features;  we  use  4149.0  cm"1  as 
the  origin  for  calculating  the  HDO  rotational  energies  in  the  cooperative  transitions  reported  in 
Table  6.  The  small  peak  at  4153.25  cm'1  remains  unassigned,  but  it  is  definitely  distinguishable 
from  the  0H2  induced  Qi(0)  line  which  peaks  at  4153. 10  cm'1  (2, 7). 

V.  SUMMARY 

Isolated  H2O,  D2O  and  HDO  molecules  exist  as  nearly  free  rotors  in  solid  pH2,  with 
rotational  constants  typically  within  a  few  percent  of  the  gas  phase  values.  The  linewidths  of 
H2O  transitions  indicate  that  rovibrational  dephasing  of  the  upper  states  in  Type  B  bands  by 
dopant-host  interactions  is  much  faster  than  dephasing  of  the  upper  states  in  Type  A  bands  by 
IVR.  Lineshape  changes  upon  annealing  indicate  that  inhomogeneous  broadening  of  the  water 
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monomer  transitions  is  of  order  0.1  cm'1  in  RVD  pHb  solids.  The  gas-to-matrix  red  shifts  of  the 
water  monomer  rovibrational  transitions  increase  monotonically  with  increasing  excitation 
energy,  but  the  details  remain  poorly  understood.  NSC  of  isolated  0H2O  and  PD2O  monomers 
follow  first  order  kinetics,  with  single  exponential  decay  lifetimes  at  T  =  2.4  K  of  1900±100  s, 
and  860±50  s,  respectively. 

Formation  of  oH2-water  complexes,  even  at  residual  0H2  concentrations  of  only  100  ppm, 
results  in  the  appearance  of  a  large  number  of  additional  absorption  features.  The  main  features 
of  the  0H2-H2O  and  0H2-D2O  spectra  are  explained  qualitatively  by  assuming  a  semi-rigid  C2V 
structure  with  the  0H2  acting  as  a  proton  donor  to  the  O  atom.  NSC  of  oFE-water  complexes 
proceeds  at  very  nearly  the  same  rate  as  for  the  corresponding  water  monomer,  refuting  the 
proposed  role  of  pairwise  magnetic  dipole-dipole  interactions  between  0H2  and  J  =  1  water 
molecules  in  accelerating  the  NSC  process.  At  elevated  0H2  concentrations  of  1 000  ppm, 
absorptions  of  larger  (oH2)n-water  clusters  are  observed.  Somehow  0H2O  and  PD2O  molecules 
clustered  with  several  0H2  molecules  survive  longer  than  the  isolated  J  =  1  water  monomers. 

The  presence  of  water  dopants  results  in  a  number  of  dopant-induced  IR  absorption  features 
of  the  pH2  host,  and  in  cooperative  water-pFb  transitions  in  which  the  vibrational  excitation  of 
the  pH2  solid  is  accompanied  by  a  pure  rotational  transition  of  the  water  dopant.  The  linewidths 
of  these  pure  rotational  transitions  match  well  with  those  for  the  direct  IR  Type  B  rovibrational 
transitions  of  the  water  dopants,  indicating  that  rotational  dephasing  is  the  dominant 
linebroadening  mechanism  in  both  cases. 
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Table  1.  Absorptions  of  H20  monomers  and  oH2-H20  complexes  at  T  =  2.4  K  in  pH2  solids  containing  ~  100  ppm 
oH2.  Peak  positions  (FWHM)  and  energy  differences  given  in  cm'1. 


■as-deposited - 1  | - annealed- 


assignment 

gas  Dhasea 

in  dH- 

AEfgas-oHA 

in  dH? 

AEfgas-DHA 

oH2-H20 

1593.7  mult. 

1593.7  mult. 

oH2-oH20 

1606.78 

n.o. 

oH2-pH20 

1608.50 

1608.48 

v2  [1  io  loi] 

1616.712 

1612.0  (7.4) 

4.7 

n.o. 

< 

to 

r— i 

t 

o 

o 

p 

1634.967 

1631.6(1.7) 

3.4 

1631.5(1.6) 

3.5 

v2  [2i2  <—  loi] 

1653.267 

1649.0  (3.0) 

4.3 

n.o. 

2v2  [In  —  O0o] 

3196.093 

3190.7  w 

5.4 

3190.7  (1.5) 

5.4 

oH2-H20 

3645.4  mult. 

3645.4  mult. 

V1  [1 10  loi] 

3674.697 

3660.5  (1.4) 

14.2 

n.o. 

vi  [lu  Ooo] 

3693.294 

3681.2(2.7) 

12.1 

3681.2  (2.2) 

12.1 

< 

73 

io 

t 

o 

U—l 

3711.103 

3698.4  (3.1) 

12.7 

n.o. 

v3  [Ooo  <—  loi] 

3732.135 

3719.79  (0.24) 

12.35 

n.o. 

oH2-oH20  [0,0  <-  1,1] 

3725.46 

n.o. 

oH2-oH20 

3728.33 

n.o. 

oH2-PH20  [1,1  <-  0,0] 

3758.46 

3758.44 

1 

v3  [loi  < —  Ooo] 

3779.493 

3765.50  (0.23) 

13.99 

3765.47(0.18) 

14.02 

oH2-pH20 

3773.44 

3773.42 

V3  P02  <-  loi] 

3801.420 

3787.1  (2.5) 

14.32 

n.o. 

V2+V3  [Ooo  < —  loi] 

5307.475 

5292.83  (0.26) 

14.65 

n.o. 

oH2-oHzO  [0,0  <-  1,1] 

5299.2  w 

n.o. 

oH2-pH20  [1,1  <- 0,0] 

5332.3  w 

5332.26 

v2+v3  [l0i  Ooo] 

5354.873 

5338.58  (0.23) 

16.29 

5338.54  (0.14) 

16.33 

v2+v3  [202  <—  loi] 

5376.945 

5360.3  (2.0) 

16.6 

n.o. 

2v2+v3  [Ooo  loi] 

6847.726 

6831.20(0.23) 

16.53 

n.o. 

2v2+V3  [l0l  Ooo] 

6895.148 

6876.91  (0.27) 

18.24 

6876.91  (0.20) 

18.24 

V1+V3  [00o  loi] 

7226.024 

7203.4(0.5) 

22.6 

n.o. 

Vl+V3  [loi  <-  Ooo] 

7273.000 

7248.69  (0.23) 

24.31 

7248.65  (0.12) 

24.35 

V1+V3  [2q2  loi] 

7294.136 

7269.4  (2.2) 

24.7 

n.o. 

a  =  Refs.  (58,59) 
mult.  =  multiple  peaks 
n.o.  =  not  observed 
w  =  weak 
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Table  2.  Absorptions  of  D20  monomers  and  oH2-D20  complexes  at  T  =  2.4  K  in  pH2  solids  containing  ~  100  ppm 
oH2.  Peak  positions  (FWHM)  and  energy  differences  given  in  cm'1 . 


■as-deposited - 1  | - annealed- 


assignment 

gas  Dhasea 

in  nHz 

AE(gas-pH,) 

in  dH? 

AE(gas-oH-)) 

1 177.4  mult. 

1177.4  mult. 

V2  [1,0  lo,] 

1190.221 

1188.10(0.27) 

2.12 

1188.07  (0.16) 

2.15 

V2  [111  Ooo] 

1199.793 

1197.71  (0.7) 

2.08 

1197.70  (0.7) 

2.09 

v2  [2,2  <—  loi] 

1209.388 

1207.1  (1.5) 

2.3 

n.o. 

2v2  [In  <—  Ooo] 

2359.689 

2356.4  w 

3.3 

2356.4  w 

3.3 

2663.4  mult. 

2663.4  mult. 

V1  [1 10  <“  loi] 

2681.894 

2673.35  (0.22) 

8.54 

n.o. 

vi  [111  Ooo] 

2691.606 

2683.0(1.1) 

8.6 

2682.9(1.1) 

8.7 

v,  [2,2  4-  loi] 

2701.004 

2692  w 

9 

n.o. 

V3  [Ooo  loi] 

2775.602 

2765.93  (0.25) 

9.67 

2765.97  (0.10) 

9.63 

oH2-pD20  [0,0  <-  1,1] 

2769.27 

n.o. 

oH2-pD20 

2771.2  w 

n.o. 

oH2-oD20 

2782.62 

2782.59 

oH2-oD20  [1,1  <-  0,0] 

2785.22 

2785.05' 

2785.15' 

2785.29' 

v3  [loi  Ooo] 

2799.759 

2789.43  (0.22) 

10.33 

2789.39  (0.11) 

10.37 

oH2-D20 

2793.8  w 

2793.8  w 

V3  [20Z  <—  loi] 

2811.214 

2800.4(1.6) 

10.8 

n.o. 

V2+V3  [00o  <—  loi] 

3943.896 

3932.55  (0.38) 

11.35 

n.o. 

oH2-oD20 

3952.25 

3952.10' 

3952.20' 

3952.33' 

V2+V3  [loi  4—  Ooo] 

3968.069 

3956.07  (0.24) 

12.00 

3956.02(0.12) 

12.05 

V2+V3  [202  4—  loi] 

3979.562 

3966.8  w 

12.8 

n.o. 

V1+V3  [loi  4-  Ooo] 

5385.788 

5367.73  (0.24) 

18.06 

5367.67(0.14) 

18.12 

a  =  Refs.  (57,61,62) 
mult.  =  multiple  peaks 
n.o.  =  not  observed 
t  =  components  of  a  triplet 
w  =  weak 
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Table  3.  Absorptions  of  HDO  monomers  and  oH2-HDO  complexes  at  T  =  2.4  K  in  pH2  solids  containing  ~  100  ppm 
oH2.  Peak  positions  (FWHM)  and  energy  differences  given  in  cm'1. 


■as-deposited - 1  | - annealed- 


assignment 

gas  nhasea 

in  nH-> 

AE(aas-DH?) 

innfL  AEfaas-nHA 

v2  [00o  loi] 

1387.976 

n.o. 

1386.91h  (0.12) 

1.07 

oH2-HDO 

1402.5  mult. 

1402.5  mult. 

oH2-HDO 

1413.28 

1413.20d 

1413.29d 

v2  [loi  <—  Ooo] 

1419.041 

1416.96  (0.27) 

2.08 

1416.92(0.17) 

2.12 

V2  [111  Ooo] 

1435.300 

1432.5(1.8) 

2.8 

1432.6(1.5) 

2.7 

V2  [2j2  <—  loi] 

1447.997 

n.o. 

1445 ,0h  (1.4) 

3.0 

vi  [Ooo  loi] 

2708.172 

n.o. 

2700.25h  (0.08) 

7.92 

oH2-HDO 

2714.90 

2714.90 

oH2-HDO 

2725.5  mult. 

2725.48' 

2725.51' 

2725.58' 

o 

o 

O 

i 

o 

> 

2738.914 

2730.02  (0.15) 

8.89 

2729.995d  (0.03) 

8.919 

2730.085d  (0.02)  w 

8.829 

2v2  [loi  Ooo] 

2797.610 

2793.82  (0.26) 

3.79 

2793.79(0.18) 

3.82 

2v2  [1 1 1  < —  0oo] 

2816.214 

2811.6  (1.2)  w 

4.6 

2811.6  (1.2)  w 

4.6 

v3  [Ooo  loi] 

3691.958 

n.o. 

3680. 37h  (0.07)  w 

11.59 

oH2-HDO 

3702.55 

3702.55 

oH2-HDO 

3705.5  mult. 

3705.39' 

3705.44' 

3705.53' 

v3  [loi  Ooo] 

3722.886 

3710.27  (0.23) 

12.62 

3710.23  (0.11) 

12.66 

V3  [In  < —  Ooo] 

3736.160 

3723.1  (1.4) 

13.1 

3723.1  (1.4) 

13.1 

V2+V3  [loi  Ooo] 

5105.027 

5090.4  w 

14.6 

a  =  Refs.  (59,60) 
d  =  component  of  doublet 
h  =  hot  band  (T  =  4.8  K) 
mult.  =  multiple  peaks 
n.o.  =  not  observed 
t  =  component  of  triplet 
w  =  weak 
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Table  4.  Gas-phase  vs.  solid-pH2  excited  vibrational  state  rotational  spacings  and  rotational  constants  sums.  Energy 
differences  and  rotational  constants  given  in  cm'1. 


vibrational 

mode 

rotational 
level  splitting 

rotational 
constants  sum 

gas  phase3 

solid  pH2 

pH2/gas 

v2  H20 

K> 

to 

1 

0 

4C 

36.555 

37.0 

1.01 

v,  H20 

2i2  - 1  xo 

4C 

36.406 

37.9 

1.04 

V3  h2o 

2o2  •  Ooo 

4B  +  A  +  C 

69.285 

67.3 

0.97 

V2+V3  H20 

2o2  *  Ooo 

4B+A  +  C 

69.470 

67.5 

0.97 

V1+V3  H20 

2o2  -  Ooo 

4B+A  +  C 

68.112 

66.0 

0.97 

v2  D20 

2i2  -  1 10 

4C 

19.167 

19.0 

0.99 

Vi  D2O 

K> 

to 

1 

t— » 

0 

4C 

19.110 

18.7 

0.98 

V3  D2O 

2o2  ■  Ooo 

4B+A  +  C 

35.612 

34.5 

0.97 

V2+V3  D2O 

to 

0 

to 

1 

O 

0 

0 

4B  +  A  +  C 

35.666 

34.2 

0.96 

v2  HDO 

2i2  -  Ooo 

4C  +  A  +  B 

60.021 

58.1 

0.97 

v2  HDO 

In  -  loi 

A-B 

16.259 

15.5 

0.95 

2v2  HDO 

In  -  loi 

A-B 

18.604 

17.8 

0.96 

v3  HDO 

In  -  loi 

A-B 

13.274 

12.8 

0.96 

a  =  Refs.  (57,61,62) 


Table  5.  Gas-phase  vs.  solid-pH2  vibrationally  averaged  rotational  constants;  AE^m  defined  in  Eq.  [4]. 


transition 

AE^Cgas)  (cm'1)3 

AEasym(pH2)  (cm'1) 

AE(pH2)/AE(gas) 

V3  H2O 

47.358 

45.71 

0.965 

v2+v3  H20 

47.398 

45.75 

0.965 

2V2+V3  H2O 

47.422 

45.71 

0.964 

V3  D2O 

24.157 

23.50 

0.973 

V2+V3  D2O 

24.173 

23.52 

0.973 

vi  HDO 

30.742 

29.79 

0.969 

v3  HDO 

30.928 

29.86 

0.966 

a  =  Refs,  (gas  1,5, 6) 
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Table  6.  Dopant-induced  and  cooperative  absorptions  in  as-deposited  pH2  solids  at  T  =  2.4  K.  Peak  positions 
(FWHM)  and  energy  differences  given  in  cm-1. 


assignment 

in  dH7 

AEfrotl 

gas  Dhasea 

AF/Reas 

OH2-OH20 

4139.41  (0.03) 

oH2-pH20 

4139.65  (0.05) 

Qi(0)H2{pH2O} 

4149.03  (0.03) 

Qi(0)  H2  {oH20} 

4149.67  (0.1) 

Q1(0)H2{oH2O} 

4149.81  (0.1) 

+  [1 10  loi]  oH20 

4164.9(1.3) 

15.2 

18.578 

0.818 

+  [In  <—  Ooo]  pH20 

4185.2  (2.4) 

36.2 

37.137 

0.975 

+  [2U  <—  loi]  oH20 

4203  (3) 

53 

55.702 

0.95 

+  [202  <—  O0o]  pH20 

4217.5  (2.5) 

68.5 

70.091 

0.977 

oH2“D20 

4139.13  (0.08) 

oH2-D20 

4139.35  (0.07) 

oH2-D20 

4139.60  (0.07) 

Qi(0)  H2  {oD20} 

4148.97  (0.03) 

Q,(0)H2{pD2O} 

4149.61  (0.1) 

Q,(0)H2{pD2O} 

4149.76  (0.2) 

? 

4153.25  (0.2) 

+  [1 10  loi]  pD20 

4159.35  (0.2) 

9.7 

10.567 

0.92 

+  [lu  4—  Ooo]  oD20 

4169.03  (0.9) 

20.06 

20.259 

0.990 

+  Pi2  loi]  pD20 

4178.2  (1.7) 

28.5 

29.952 

0.95 

+  [2q2  <—  Ooo]  oD20 

4184.1  (1.5) 

35.1 

35.878 

0.98 

+  [loi  Ooo]  HDO 

4164.5  (0.5) 

15.5 

15.508 

1.00 

+  [I11  Oqo]  HDO 

4178.4  (1) 

29.4 

29.809 

0.99 

a  =  Refs.  (57,61,62) 
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FIGURE  CAPTIONS 


FIG.  1 .  Energy  levels  and  allowed  rovibrational  transitions  for  H2O  and  D2O  molecules;  adapted 
from  Refs.  (37,54).  Type  A  bands  have  Av3  odd,  e.g.:  v3,  V2+V3,  2v2+V3,  V1+V3.  Type  B  bands 
have  AV3  even,  e.g. :  Vi,  V2. 

FIG.  2.  IR  absorption  spectra  (V2  H2O  region)  showing  the  effects  of  0H2O  — >  PH2O  conversion. 
The  0.7-mm-thick  30  ppm  FEO/pFE  sample  was  deposited  at  T  =  2.4  K  over  a  15  minute  period. 
Trace  (a)  is  for  the  1 1-minute-old  sample  at  T  =  2.4  K;  trace  (b)  is  for  the  42-minute-old  sample 
at  T  =  2.4  K;  trace  (c)  is  for  the  71 -minute-old  sample  at  T  =  2.4  K;  trace  (d)  is  for  the  272- 
minute-old  sample  at  T  =  2.4  K,  after  annealing  at  T  =  4.8  K  for  10  minutes.  Each  spectrum 
includes  20  co-added  scans  acquired  over  a  2.0  minute  period  (centered  at  the  times  given 
above),  and  is  presented  at  0.1  cm'1  resolution.  Features  marked  “0”  are  stronger  in  samples 
deliberately  doped  with  higher  0H2  concentrations,  features  marked  are  coincident  with 
atmospheric  water  vapor  absorptions. 

FIG.  3.  IR  absorption  spectra  (V3  H2O  region)  showing  the  effects  of  oH20  ->  pH20  conversion 
for  the  sample  described  in  FIG.  2.  Trace  (a)  is  for  the  6-minute-old  sample  at  T  =  2.4  K;  trace 
(b)  is  for  the  36-minute-old  sample  at  T  =  2.4  K;  trace  (c)  is  for  the  66-minute-old  sample  at  T  = 
2.4  K;  trace  (d)  is  for  the  260-minute-old  sample  at  T  =  2.4  K,  after  annealing  at  T  =  4.8  K  for  10 
minutes.  Each  spectrum  includes  10  co-added  scans  acquired  over  a  2.0  minute  period  (centered 
at  the  times  given  above),  and  is  presented  at  0.05  cm'1  resolution.  Features  marked  “o”  are 
stronger  in  samples  deliberately  doped  with  higher  0H2  concentrations,  features  marked  are 
coincident  with  atmospheric  water  vapor  absorptions,  the  feature  marked  “i”  is  due  to  a  different 
water  isotopomer  (HDO). 

FIG.  4.  IR  absorption  spectra  (V2  D2O  region)  showing  the  effects  of  PD2O  -»  0D2O  conversion. 
The  0.6-mm-thick  30  ppm  D2O/PH2  sample  was  deposited  at  T  =  2.4  K  over  a  10  minute  period. 
Trace  (a)  is  for  the  1 -minute-old  sample  at  T  =  2.4  K;  trace  (b)  is  for  the  16-minute-old  sample  at 
T  =  2.4  K;  trace  (c)  is  for  the  31-minute-old  sample  at  T  =  2.4  K;  trace  (d)  is  for  the  185-minute- 
old  sample  at  T  =  2.4  K,  after  annealing  at  T  =  4.8  K  for  20  minutes.  Each  spectrum  includes  20 
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co-added  scans  acquired  over  a  2.0  minute  period  (centered  at  the  times  given  above),  and  is 
presented  at  0.1  cm'1  resolution.  Features  marked  “o”  are  stronger  in  samples  deliberately  doped 
with  higher  0H2  concentrations. 

FIG.  5.  IR  absorption  spectra  (V3  D2O  region)  showing  the  effects  of  PD2O  ->  oD20  conversion 
for  the  sample  described  in  FIG.  4.  Trace  (a)  is  for  the  6-minute-old  sample  at  T  =  2.4  K;  trace 
(b)  is  for  the  36-minute-old  sample  at  T  =  2.4  K;  trace  (c)  is  for  the  192 -minute-old  sample  at  T  = 
2.4  K,  after  annealing  at  T  =  4.8  K  for  20  minutes.  Each  spectrum  includes  10  co-added  scans 
acquired  over  a  2.0  minute  period  (centered  at  the  times  given  above),  and  is  presented  at  0.05 
cm'1  resolution.  Features  marked  “0”  are  stronger  in  samples  deliberately  doped  with  higher  0H2 
concentrations. 

FIG.  6.  Kinetics  of  0H2O  ->  pFEO  and  pD2<D  ->  0D2O  conversion  at  T  =  2.4  K.  The  solid  lines 
are  single  exponential  fits  to  the  decay  of  the  J=1  water  species. 

FIG.  7.  IR  absorption  spectra  (V2  HDO  region)  showing  the  effects  of  temperature  cycling.  The 
1.8-mm-thick  H2O/HDO/D2O/PH2  sample  containing  5  ppm  HDO  was  deposited  at  T  =  2.4  K. 

In  sequence:  trace  (a)  T  =  2.4  K;  trace  (b)  T  =  4.8  K;  trace  (c)  T  =  2.4  K;  trace  (d)  T  =  4.8  K; 
trace  (e)  T  =  2.4  K.  Each  spectrum  is  presented  at  0.01  cm'1  resolution.  Features  marked  “0”  are 
stronger  in  samples  deliberately  doped  with  higher  0H2  concentrations. 

FIG.  8.  IR  absorption  spectra  (V3  HDO  region)  showing  the  effects  of  temperature  cycling  for 
the  sample  described  in  FIG.  7.  In  sequence:  trace  (a)  T  =  2.4  K;  trace  (b)  T  =  4.8  K;  trace  (c)  T 
=  2.4  K;  trace  (d)  T  =  4.8  K;  trace  (e)  T  =  2.4  K.  Each  spectrum  is  presented  at  0.01  cm'1 
resolution.  Features  marked  “0”  are  stronger  in  samples  deliberately  doped  with  higher  0H2 
concentrations;  the  features  marked  “i”  are  due  to  a  different  water  isotopomer  (H2O). 

FIG.  9.  IR  absorption  spectra  (vj  HDO  region  detail)  showing  the  effects  of  temperature  cycling. 
The  1.5-mm-thick  D20/HD0/pH2  sample  containing  1  ppm  HDO  was  deposited  at  T  =  2.4  K.  In 
sequence:  trace  (a)  T  =  2.4  K;  trace  (b)  T  =  4.8  K;  trace  (c)  T  =  2.4  K.  Each  spectrum  is 
presented  at  0.01  cm'1  resolution. 
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FIG.  10.  Gas-to-matrix  peak  shifts  for  rovibrational  transitions  of  water  monomers  in  solid  pHa 
at  T=2.4  K  Data  for  as-deposited  peak  positions  listed  in  Tables  1, 2  and  3.  Points  for 
transitions  involving  V2  excitation  are  circled. 

FIG.  1 1.  ER  absorption  spectra  (V2  H2O  origin)  showing  the  effects  of  0H2/H2O  clustering. 

Trace  (a)  is  for  an  as-deposited  1.0-mm-thick  10  ppm  H2O/PH2  sample  containing  a  100  ppm 
0H2  at  T  =  2.4  K.  Trace  (b)  is  for  an  as-deposited  3.3-mm-thick  10  ppm  H2O/PH2  sample 
containing  a  1000  ppm  0H2  at  T  =  2.4  K.  Each  spectrum  is  presented  at  0.01  cm'1  resolution. 

FIG.  12.  IR  absorption  spectra  (V3  H2O  region)  showing  the  effects  of  0H2/H2O  clustering. 

Trace  (a)  is  the  same  data  shown  above  in  FIG.  3(a);  this  sample  contains  a  100  ppm  0H2.  Trace 
(b)  is  for  the  sample  described  above  in  FIG.  1 1(b);  this  sample  contains  a  1000  ppm  oH2.  Each 
spectrum  is  presented  at  0. 1  cm'1  resolution. 

FIG.  13.  IR  absorption  spectra  (induced  IR  region)  showing  the  effects  of  oH20  ->  pHtO 
conversion  for  the  sample  described  in  FIG.  2.  Traces  (a),  (b),  and  (c)  are  from  the  same  spectra 
presented  in  FIG.  3  traces  (a),  (b),  and  (c).  The  labels  “Qi(0)  H2  {X}”  are  for  pure  Av  =  +1 
vibrational  transitions  of  pHh  induced  by  species  X.  The  labels  “+  [J’Ka’Kc’  <-  J’’Ka”Kc”]  Y”  are 
shorthand  for  the  Qi(0)  H2  +  [J’Ka’Kc’  <-  J’’Ka”Kc”]  Y  cooperative  absorptions  in  which  the 
neighboring  dopant  species  Y  undergoes  a  pure  rotational  transition. 

FIG.  14.  IR  absorption  spectra  (induced  IR  region)  showing  the  effects  of  PD2O  -»  0D2O 
conversion  for  the  sample  described  in  FIG.  4.  Traces  (a)  and  (b)  are  from  the  same  spectra 
presented  in  FIG.  5  traces  (a)  and  (b).  The  spectroscopic  labels  are  explained  in  FIG.  13. 
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FIG.  1.  Fajardo,  Tam,  andDeRose. 
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FIG.  4.  Fajardo,  Tam,  and  DeRose. 
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FIG.  6.  Fajardo,  Tam,  and  DeRose. 


Approved  for  public  release;  distribution  unlimited. 


log-iodc/D 


1380  1390  1400  1410  1420  1430  1440  1450 

wavenumber  (cm'1) 


FIG.  7.  Fajardo,  Tam,  and  DeRose. 
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FIG.  9.  Fajardo,  Tam,  and  DeRose. 
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FIG.  10.  Fajardo,  Tam,  and  DeRose. 
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